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SECTION  ! 

INTRODUCTION  AND  SUMMARY 


A  localized  cloud  containing  gases  that  have  infrared  absorption  bands, 
e.g.,  the  chemical  effluent  of  a  stationary  source,  ran  be  detected  on  the 
basis  of  the  contrast  it  produces  in  an  IR  scene  image.  Obvious  requirements 
in  this  method  of  remote  detection  are  that  the  cloud  temperature  is  different 
from  the  background  brightness  temperature  and  that  *he  resultant  spatial 
radiance  contrast  is  distinguishable  from  background  clutter. 

Resolution  of  the  IP.  image  into  many  spectral  elements  will  enhance  the 
ability  of  the  system  to  distinguish  between  different  target  species  having 
overlapping  spectral  bands  and/or  to  suppress  the  effects  of  interferences, 
particu'arly  those  resulting  from  spatial  and  temporal  variations  in  atmos¬ 
pheric  spectral  absorption  and  background  emission. 

The  detection  concept  addressed  by  this  study  is  pictured  in  Fig.  1.  A 
scene  that  includes  a  chemical  effluent  cloud  is  viewed  by  an  airborne  .^pee- 
trometer  having  a  mosaic  or  detectors,  some  of  which  see  the  target  cloud- 
plus-background  and  others  that  see  only  the  background.  The  difference  in 
the  outputs  of  iwo  appropriately  selected  detectors  gives  the  spectral 
radiance  contrast  between  the  target  cloud  and  the  background.  Comparison  ot 
the  measured  contrast  spectrum  and  a  computed  reference  contrast  spectrum  can 
be  used  to  identify  the  presence  of  a  particular  trace  gas,  and  also  to 


Figure  1.  Detection  scenario 
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estimate  its  "detectable  quantity".  Detection  and  quantification  involve 
s i nip  1 e  1  east -squares  estimation  procedures  developed  in  a  previous  study' 
and  i  i-v  levied  In  '.tic.  lion  it.  The  detectable  (observable)  quantity  erf  the 
taryet  trace  gas  is  essentially  the  product  of  its  molecular  column  thickness 
in  the  i i ne-of-s i ght  and  the  Planck  spectral  radiance  difference  corresponding 
to  the  temperature  difference  between  the  target  gas  and  background.  The 
minimum  detectable  quantity  (MDQ)  will  be  defined  as  the  value  that  corres¬ 
ponds  to  approximately  95  percent  probability  of  detection  and  one  percent 
probability  of  false  detection.  The  MDQ.  depends  on  the  noise  equivalent 
spectral  radiance  (NESR)  of  the  spectrometer  system,  the  variance  of  the 
reference  contrast  spectrum  and  the  number  of  resolved  spectral  elements  in 
the  detection  band. 

Of  course,  in  some  applications  it  may  be  desirable  or  necessary  to 
remotely  infer  both  the  target  gas  column  thickness  and  its  temperature.  In 
principle,  these  quantities  can  be  determined  by  more  detailed  analysis  of 
the  observed  contrast  spectral  radiance,  but  it  is  not  possible  to  specify 
detection  limits  and  accuracies  without  further  study  of  the  problem. 

Knowledge  of  the  predicted  MDQ  can  be  used  at  least  to  determine  whether 
particular  combinations  of  target  gas  column  thickness  and  gas-background 
temperature  difference  will  be  detectable  in  a  given  scenario. 

A  previous  study'  by  Atmospheric  Radiation  Consultants,  Inc.  (ARC) 

established  baseline  design  configurations  for  an  airborne  Fourier  Transform 

Spectrometer  (FTS)  detection  system  with  modest -s i zed  foreoptics.  The  study 

determined  the  NESR  of  one  of  the  baseline  systems  and  the  corresponding 

MDQ's  for  the  15  species  C^C^,  CH^,  CH^I,  CO,  C02>  DF,  HBr,  HC1,  HF ,  HI, 

HNO^,  NH^,  N2O.  NO^  and  SO^.  The  MDQ  results  given  in  Ref.  1  assume  that 

the  target  cloud  is  viewed  in  the  direction  of  the  nadir  from  an  altitude 

above  most  of  the  absorbing  atmosphere.  (These  results  and  the  detection 

'>  ? 

concept  are  also  summarized  in  two  papers;  cf.  Zachor,  ct  a  !*' ' J)  . 

The  choice  of  an  FTS  sy.tem  over  a  grating  spectrometer  system  for  this 
application  is  supported  by  a  study  performed  by  SSG,  Inc.'4  The  basic 

radiometric  design  requirements  and  tradeoffs  were  established  in  an  earlier 

r 

study""  by  Bartlett  Systems,  .nc, 
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The  PROOF  OF  CONCEPT  TESTING  program  consisted  of  1 aboratory/f  i e I d  tests 
and  data  analyses  by  Utah  State  University,  ARC  and  the  Aerospace  forpoiation. 
The  purpose  of  the  USU  Electrodynamics  laboratory  and  ARC  efforts  was  to 
validate  the  detection  concept  and  theoretical  prediction  models  of  Ref.  1 
through  field  tests  and  subsequent  data  analysis,  as  we i 1  as  through  computer 
simulation.  ARC  was  also  charged  with  investigating  methods  for  obtaining 
separate  estimates  of  the  target  gas  column  thickness  and  temperature. 


The  field  tests  were  performed  by  USU,  with  guidance  from  AFGL 
personnel.  A  model  stack  emitting  a  controlled  flow  of  ambient-temperature 
N^O  in  front  of  a  heated  "background"  plate  was  observed  by  an  FTS  system  at 
a  distance  of  273  meters.  The  source  and  detection  system  were  both  at  ground 
level.  The  side  walls  of  the  "three-wall"  model  stack  limited  the  expansion 
of  the  stack  effluent  and  reduced  the  effects  of  wind.  It  was  determined  by 
USU  that  the  N20  concentration  {and  column  thickness)  observed  by  the  FTS 
system  could  be  reliably  controlled  by  metering  the  N^O  supply  to  the  stack. 

1 nterferograms  were  recorded  for  M^O-pl us-background  (supply  ON)  and  for 
background  only  (supply  OFF).  USU  reduced  the  i nf erferograms  to  spectra  and 
supplied  these  to  ARC  for  analysis.  They  also  supplied  appropriate  meteorolo¬ 
gical  data  and  the  background  plate  temperature.  From  these  data  ARC  computed 
the  reference  contrast  spectrum,  the  system  NESR,  the  MDQ  and  corresponding 
detection  threshold  quantity,  the  detection  decision  (vos  or  no),  the  detect¬ 
able  quantity  present  and  the  corresponding  volumn  concentration  of  N„,C. 
Afterwards,  the  inferred  concentrations  were  compared  to  values  supplied  by  USD. 


The  detection  band  selected  for  the  experiment  is  21/0-2260  cm 


1 


wh  i  ch 


includes  most  of  the  band  of  f^O  (this  is  a  good  choice  for  the  field  tests, 
but  not  the  optimum  oectral  band  for  detection  of  N^O  from  a  high-altitude 


J  -J  |  £  4 

platform').  Atmospheric  N70  and  the  isotope  C  0^  produce  significant 

t-  *-• 

absorptions  over  a  273-meter  ground- level  path  .n  the  selected  band.  Thus, 
the  tests  included  the  effects  of  atmospheric  attenuation,  even  though  the 
path  length  was  short.  Detection  near  the  MDQ  level  implies  very  low  signal 
to  noise  in  the  measured  contrast  radiance  spectrum.  The  field  tests  were 
designed  primarily  to  verify  that  reliable  detection  could  be  achieved  at,  or 
slightly  above,  the  MDQ  level,  in  a  benign  real-world  environment. 


The  computer  simulation  performed  by  ARC  utilized  a  simile  pair  of  IISU 
measurements  representing  a  contrast  radiance  spectrum  of  high  s i yna I - to-no i so . 

I  1 1  n  i  t  mi  1 111 1  «- 1  wn  u  s  r  1 1  ri-.rnl  I  nl  I  y  t  <  i  mid  psr  i  till  i  ■  r  nfwli  mi  rn  >  I  s  r  ,  I  .  e  .  ,  to  dr 1 1  r  a  dr 

the  S/N,  or  equivalently,  to  raise  the  MOO.  level  to  the  detectable  quantity 
present  during  the  HSU  high-S/N  measurement.  The  resultant  noisy  contrast 
spectra  were  analyzed  in  the  same  way  as  contrast  spectra  obtained  from  pairs 
of  USU  low-S/N  measurements.  The  noi se~add i t i on  process  was  performed  on  the 
computer  one  thousand  times.  From  the  number  of  successful  detections  it  was 
verified  that  the  detection  probability  was  approximately  95  percent.  Similar 
computations  involving  noise-only  spectra  were  used  to  verify  that  the  false 
detection  probability  was  approximately  one  percent. 

This  final  report  summarizes  ARC'S  efforts  in  the  PROOF  OF  CONCEPT  TESTING 
program.  Details  of  the  field  tests,  the  three-wall  stack,  the  ITS  system 
used  in  the  tests,  reduction  of  i nferferograms ,  etc.  are  contained  in  USU's 
final  report. 

The  remainder  of  Section  1  gives  a  sample  of  the  results  obtained  and  the 
conclusions  of  the  ARC  study.  Section  2  reviews  the  theoretical  basis  of  the 
concept  and  the  spectral  analysis  methods.  Sections  3  through  5  contain  addi¬ 
tional  results  and  a  more  detailed  account  of  the  simulation  and  data  analysis 
efforts . 

1 . 1  Sample  of  Results 

The  spectral  radiance  data  for  each  field  test  was  supplied  by  USU  as  a 
digital  tape  which  usually  contained  12  consecutive  target  spectra  (N^O  supply 
on),  12  corresponding  background  spectra  (N^O  supply  off),  and  two  spectral 
responsivity  runs.  The  N^O  flow  rate  was  held  constant  in  each  series  of  12 
target  runs,  except  during  the  test  of  19  February  1981,  when  it  was  purposely 
varied. 

Figure  2  compares  the  N^O  concentrations  that  were  deduced  by  ARC  from 
spectral  analysis  of  the  19  February  data  to  those  measured  by  USU  (from  the 
calibrated  flow  meter).  The  analysis  procedure  is  capable  ol  providing  esti¬ 
mates  of  only  the  detectable  quantity,  but  knowledge  ot  the  ambient  ( N 2 0 ) 
temperature  and  pressure,  the  background  temperature  and  physical  length  of 


k 


the  I i ne- of - s i ght  through  the  N.,0  allowed  n  conversion  to  N,0  volunm  conc.cn- 

<t  i. 

tration  (ppmV  units).  It  is  seen  that  the  actual  concentration  decreased 

between  runs  I  and  6,  then  underwent  a  symmetrical  increase  between  runs  5  and 

11,  and  finally  decreased  from  run  11  to  run  17.  The  figure  shows  that  N.,0 

was  successfully  detected  in  the  eight  runs  for  which  the  actua1  concentration 

was  above  the  theoretical  detection  threshold;  in  the  remaining  four  cases  it 

was  not  detected  (there  were  no  false  detections).  Moreover,  the  differences 

between  the  deduced  and  measured  values  are  within  the  theoretical  '  one-sigma 

uncertainty  in  the  deductions.  Note  that  two  successful  detections  were  made 

10  2  2  - 1 

very  near  the  MDQ  level  of  17.8  ppi nV  =  1.4  x  10  (mol ec/cm  ) (W/cm  sr  cm  ). 

The  constant  flow  rates  used  by  USU  in  the  other  series  of  measurements 
corresponded  to  concentrations  ranging  from  approx i mat e 1 y  1.5  to  4  times  the 
MDQ  concent  rat i on .  Very  good  agreemen  between  inferred  and  measured  concen¬ 
trations  was  also  obtained  f or  these  cases. 


1  he  estimate  of  theoretical  uncertainty  may  be  high,  since  it  is  proportional 
to  the  system  NESK,  which  was  estimated  from  the  difference  between  successive 
background  spectra  Thot  is,  temporal  background  variations  may  have  caused 
the  uncertainty  to  be  overestimated.  This  r.^y  explain  the  near -symmet ry 
(approximate  repeatability)  of  the  deduced  concentrations .  It  is  possible 
that  the  deduced  concent  rat  i  ons  j'-e  mor--  accurate,  than  the  measured  values. 
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F inure  3  shows  a  computed  reference  contrast  spectrum  for  one  of  the 
data  runs,  and  two  examples  of  measured  contrast  spectra  that  resulted  in  N^O 
detection.  The  middle  panel  of  the  future  shows  the  measured  contrast  when 
the  detectable  amount  is  approximately  four  times  the  MDQ.  At  this  level  the 
presence  of  N?0  is  barely  discernible  by  visual  examination  of  the  spectrum. 

In  tiie  bottom  panel  the  detectable  amount  is  —  1.5  times  the  MDQ.  There  is 
extremely  high  probability  of  detection  (>s  0.95)  when  the  detectable  amount 
is  1.5  times  the  MDQ,  even  though  the  measured  contrast  spectrum  for  this 
amount  looks  like  noise.  The  spectrum  s i gna 1  - to-noi se  is  only  0.26  for  the 
bottom  panel  in  Fig.  3;  at  the  MDQ  level  the  S/N  would  be  0 . 1 5 .  The  "flat 
spots"  between  2250  and  2260  cm  '  in  the  two  measured  contrast  spectra  repre¬ 
sent  spectral  elements  that  were  excluded  in  the  correlation  detection 
processing.  The  excluded  elements  have  very  low  atmospheric  transmittance 
(hence,  little  information)  due  to  the  very  strong  isotopic  CO^  lines  seen  in 
the  top  panel  of  the  figure. 

1 . 2  Cone1 us i ons 

Analysis  of  the  field  test  data  revealed  some  unanticipated  problems  both 
with  the  data  analysis  methods  and  the  FTS  instrumentation.  There  was  suffi¬ 
cient  time  remaining  in  the  program  after  these  problems  were  resolved  to 
obtain  valid  data  and  to  perform  valid  analyses  of  the  data.  The  identification 
and  solution  of  the  problems,  which  are  discussed  in  Sections  3  and  k ,  was  a 
useful  byproduct  of  the  program. 

The  major  conclusions  of  the  study  are: 

(1)  The  concept  has  been  validated  completely  by  computer  simulation. 

The  simulation  verified  the  predicted  detection  probability,  false 
detection  probability,  and  uncertainty  in  the  detected  quantity  at 
levels  equal  to  and  lower  than  the  MDQ  level.  The  simulations 
used  a  field-measured  contrast  spectrum,  and  computer-generated 

no i se . 

(2)  Detectable  amounts  obtained  by  spectral  analysis  of  the  field  data 
taken  in  January  and  February  are  in  excellent  agreement  with 


Actually  -ANV  =  1  s  shown.  ANV  is  negative  because  N^O  was  observed  in 

absorption  rather  than  emission. 
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Figure  3.  Normalized  reference  contrast  spectrum  for  vj  band  of  N?0  (top  panel) 
and  two  normalized  contrast  spectra  obtained  from  measurements. 


directly  measured  amounts.  The  quantity  of  field  data  obtained 
was  not  sufficient  to  allow  estimation  of  a:tual  detection  and 
false  detection  probabilities. 

(3)  The  calculated  MDQ's  and  thresholds  are  possibly  too  high,  as  a 
result  of  unknown  background  temperature  variations.  Tu : s  means 
only  that  we  may  have  underestimated  the  detection  capabilities 
of  the  USU  FTS  field  system;  the  possible  error  in  MDQ  does  not 
affect  the  inferred  amounts. 

(h)  We  examined  the  possibility  of  inferring  the  target  gas  column 
thickness  and  temperature  and  the  background  temperature  by  a 
more  detailed  analysis  of  the  measured  spectra,  and  have  concluded 
that  these  quantities  can  be  determined  with  reasonable  accuracy 
when  the  detectable  quantity  is  much  greater  than  the  MDQ, 

However,  if  the  target  cloud  contains  another  species  (such  as  CO^) 
that  may  not  be  of  particular  interest  but  is  present  in  suffi¬ 
ciently  larqe  quantity,  the  cloud  temperature  can  be  determined 
with  high  accuracy  by  a  detailed  analysis  of  its  spectrum,  and  then 
the  target  gas  column  thickness  can  be  inferred  with  high  accuracy, 
even  if  the  target  gas  is  present  at  the  MDQ  level.  The  capabili¬ 
ties  of  the  method  to  infer  target  gas  column  thickness  and  tempera¬ 
ture  requires  a  detailed  specification  of  the  expected  total 
composition  of  the  cloud. 


SECTION  2 
REVIEW  OF  THEORY 


The  idea  of  a  minimum  detectable  quantity  implies  low  target  gas 
concentrations.  It  also  implies  a  benign  background,  since  background  (and 
atmospheric)  interferences  would  generally  require  using  a  higher  detection 
threshold,  which  would  allow  detection  only  at  higher  concentration  levels. 
Specifically,  it  will  be  assumed  that  the  target  cloud  is  optically  thin  in 
the  spectral  bands  used  for  detection,  and  that  the  background  and  atmosphere 
in  the  two  fields  of  view  (IFOV's)  used  to  obtain  the  spectral  contrast  are 
identical.  The  only  type  of  interference  included  in  this  section  is  system 
(detector)  noise.  The  effects  of  some  types  of  atmospheric  and  background 
interferences  are  analyzed  in  Ref.  1. 


Under  the  stated  assumptions  the  spectral  radiance  contrast  for  the  two 
IFOV's  indicated  in  Fig.  1  is  given  by' 


AN  =  N. 


N ,  —  Dt  o_  ; 

1  V  cv 


D  ■ 


in 


where 

t  =  atmospheric  spectral  transmittance  between  the  target  and 
sensor , 

spectral  absorption  coefficient  of  the  (single)  target  gas, 

molecular  column  thickness  of  the  target  gas  (number  of 
molecules  per  unit  area  in  the  1 i ne-of-s i ght) , 

and  AB  =  difference  between  the  spectral  radiance  of  a  blackbody  at 
the  gas  temperature  Tg  and  the  actual  spectral  radiance  of 
the  background  at  temperature  Tsfc,  averaged  over  the  detec¬ 
tion  spectral  band. 


The  detectable  quantity,  denoted  by  D,  is  essentially  a  scale  factor  in  the 
measured  contrast  spectrum.  Equation  (1)  is  an  approximation  only  because  we 


have  replaced  the  actual  spectral  variation  AB^  by  its  average  value. 
However,  the  variation  is  slow  compared  to  that  of  and  ot  ^  and  is  nearly 
linear  over  narrow  detection  bands,  so  that  the  approximation  is  a  good  one, 
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especially  If  xy  —  1  and  a  Is  more-or-less  symmetrical  about  the  center  of 
the  detection  band. 


A  reasonably  accurate  theoretical  estimate  of  x  a  can  be  obtained 

v  gv 

using  a  line-by-line  computer  code  such  as  AFGL's  FASCODI.  Then  an  estimate 
D'  of  D  can  be  obtained  by  finding  the  value  that  minimizes  the  mean  square 


Of  course,  a  is 
gv 


difference  between  the  measured  spectrum  AN  and  D"x  a 

v  v  gv 

a  f unct ■ on  of  the  gas  temperature,  but  Dina  /3T  <<  8lnAB  /9T  .  Hence,  the 

gv  g  v  g 

u1  a  of  some  guess  temperature  in  computina  a  will  not  result  in  a  large 
error  in  D'. 


Some  simple  types  of  background  and  atmospheric  interferences  can  be  sup¬ 
pressed  by  subtractina  the  means  of  AN  and  x  a  from  these  spectra  before 

v  v  gv 

they  are  used  to  estimate  D.  The  least-squares  estimator  of  D  is  then 


£(a-S)(b-b)  oab 

St-*2  5  -2 


where  a  =  AN  +  NOISE  ,  b  =  x  a  ,  and  the  summations  are  over  the  digitized 
v  v  v  gv 

spectrum  values.  Equation  (2)  is  the  prescription  for  computing  the  estimate 
D"  given  the  measurement  a  and  reference  contrast  spectrum  b.  By  combining 
Eqs .  (1)  and  (2),  we  find  that  D'  can  also  be  expressed  by 


D'  -  D  +  (3) 

°b 

where  is  the  covariance  of  the  digitized  reference  spectrum  and  the  spec¬ 

trum  noise,  and  a ^  is  the  variance  of  the  digitized  reference  spectrum. 

The  spectrum  noise  samples  will  be  Gaussian  with  zero  mean,  and  the 
samples  will  also  be  independent  provided  the  sample  spacing  is  not  smaller 
than  the  resolution  Av  of  the  FTS  system.  It  can  be  shown'  that  if  the 
spacing  is  Av  the  second  term  of  Eq.  (3)  is 


°bn  =  *°n  =  ^  Jl  NESR 

ob2  /N  ob  /M  ob 


xo 


D" 


(M 
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where  x  is  a  unit  normal  random  variate,  M  is  the  total  number  of  resolved 

spectral  elements  in  the  spectra  a  and  b,  and  NESR  is  the  average  system 

noise-equivalent  spectral  radiance  over  the  detection  band,  Note  that  the 

effective  rms  noise  o„  equals  /2  NESR  because  AN  is  the  difference  of  two 

n  v 

measurements  with  independent  noise.  Thus,  the  rms  uncertainty  in  D'  is 

aD.  =  /2  NeTr'/^  ob)  . 

The  quantity 


qbD  _  "_S_" 

Jl  NESR  N 


(5; 


is  the  rms  signal  variation  [from  Eq .  (1)  and  the  definition  of  o^]  divided 
by  the  rms  spectral  noise;  i.e.,  it  is  the  spectral  detection  eignal-to- 
noise.  Since  D  is  this  quantity  times  /M,  it  is  evident  that  for  suffi¬ 

ciently  large  M  the  detectable  quantity  can  be  estimated  accurately  even  if 
local  spectrum  features  are  well  below  the  noise  level  (S/N  <<  1). 


D  will  be  negative  if  the  target  gas  is  observed  in  absorption  (is 
"cooler"  than  the  background)  and  positive  if  it  is  observed  in  emission. 
Hence,  detection  of  the  target  gas  can  be  based  on  the  test 


O' 


7 

> 


(6) 


where  Rp  is  some  predetermined  threshold.  Knowing  the  statistics  of  D' 

(Eqs.  3  and  M  we  can  work  out  the  statistics  for  |D'j  for  both  target 
present  (D  ^  0)  and  target  not  present  (D  *  0),  and  finally  determine  a 
threshold  Rp  that  yields  satisfactory  detection  and  false  detection 
probab i 1 i t i es  ■ 

.'quations  (3)  and  (4)  show  that  the  probability  density  function  (p.d.f,) 
of  D'  is  normal  with  mean  D  and  standard  deviation  =  /2  NESR/(/m  o^) . 

The  p.f.d.  has  zero  mean  if  the  target  is  not  present.  It  is  well  known 
that  the  p.d.f.  of  the  absolute  value  of  a  random  variate  x  is 


1  1 


g !  ( I  x  j ) 


g(x)  +  g(-x) 


(7) 


if  g(x)  denotes  the  p.d.f.  of  x.  Hence,  | D "  |  has  the  p.d.f.  indicated  by 
Fig.  4  (and  the  upper  abscissa  scale).  That  is,  the  p.d.f. 's  of  |D'j  for 
D  =  0  and  D  4-  0  are  2  x  normal  and  approximately  normal,  respectively,  pro- 
vi  ed  we  can  assume  that  the  original  p.d.f.  of  D"  for  D  /  0  has  essentially 
no  overlap  with  its  mirror  image  [g(~x)  —  0  everywhere  that  g(x)  :s  appre¬ 
ciably  different  from  zero].  This  assumption  will  always  be  valid  if  the 
resulting  p.d.f. 's  can  yield  high  system  performance,  e  g.,  detection  proba¬ 
bility  equal  to  95  percent  and  false  detection  probability  equal  to  one  percent 


("S/N"  <<  I) 

Figure  h .  The  probability  density  functions  for  | D ' |  and  |p|. 
These  are  valid  for  |p|  provided  "S/N"  <<  1. 


The  detection  probability  P (T | T )  and  false  detection  probability  P(T|0) 
are  defined  by 


P (T | T) 


P  (T  |  C) 


00 


(D  +  0)  d|D'| 


(D  «■  0)  d  |  D '  |  . 


and 


(8) 


>se  probabilities  are  the  areas  under  the  two  curves  in  Fig.  fr  at  |D‘j  >  R^. 
It  is  evident  that  P (T  |  T )  depends  on  R^,  .  and  D,  and  that  P(T|0)  depends 

on  Rp  and  only.  Standard  tables  of  the  normal  distribution  (integrals  of 

a  unit  normal  p.d.f.)  can  be  used  to  show  that  the  threshold  value 


R 


D 


2 .  58 


2  58 


/2  NESR 
/M  ob 


(9) 


results  in  approximately  95  percent  detection  probability  and  approximately 
one  percent  false  detection  probability  when  D  equals 


*4 .23 


V 


A .  23 


/2  NESR 

M  ob 


=  HDQ  . 


(10) 


This,  by  definition,  is  the  minimum  detectable  quantity.  Note  that  Rq  = 

0.610  x  HDQ.. 

The  NESR ,  M  and  all  vary  with  the  spectrometer  resolution  A v .  One 
can  show  that  a  near-optimum  value  for  Av  is  the  halfwidth  of  the  target  gas 
spectral  lines,  which  i  s  ~0 . 1  cm  '  for  targets  near  sea-lt-vel.  That  is. 

Av  —  0.1  cm  *  will  give  the  smallest  MDQfromEq.  (10).  The  single  param¬ 
eter  in  Eq.  (10)  that  depends  on  the  target  gas  and  detection  geometry  is  o^, 

the  standard  deviation  of  the  reference  spectrum  to  degraded  to 
_  I  v  gv 

0.1  cm  ‘  resol ut ion . 

2 . 1  Statistics  of  the  Correlation  Coefficient 

In  the  previous  study^  a  numerical  procedure  was  developed  for  computing 
the  statistics  of  the  spectral  correlation  coefficient  p,  defined  by 


^(a-a)  (b-b) 


X>;n2!mb-b)2 


"ab 

a  a. 

a  b 


(ID 


One  can  use  spectra  having  many  points  per  Av  --  this  would  not  appreciably 
affect  the  HDQ,  but  Eq .  (10)  would  not  be  valid  for  calculating  the  HDQ 
unless  H  is  redefined.  Equations  (M,  (9)  and  (10)  are  valid  only  if  the 
H  noise  samples  are  independent. 
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The  test  for  detection  proposed  in  Ref.  1  was 


Ip!  >  R  • 

We  will  show  that  this  test  is  equivalent  to  |D‘j  >  (Eq.  6)  provided 


S/N  <<  I  . 


Comparison  of  Eqs.  (2)  and  (11)  shows  that 


p  =  (ob/oa)D'  ■ 


Substituting  Eqs.  (3)  and  {h)  ,  which  define  D',  into  (13)  >  we  obtain 


D  ♦  jl  Is. 


where  x  is  again  a  unit  normal  random  variate.  We  recall  that  a  =  AN  + 

v 

NOISE^;  hence  p  does  not  have  a  simple  (normal)  p.d.f.,  like  D'.  However, 

in  the  case  of  very  low  S/N,  o  o  ,  and 

a  n 


p  c*  —  D' 
°n 


A  D  +  JL 


_s_"  _x_ 

N  /M 


This  result  means  that  if  S/N  is  very  low,  the  random  variate  |p[  and 
|(oj./on)D>|  have  the  same  p.d.f. ‘s,  as  indicated  by  the  two  abscissa  scales 
in  Fig.  it.  In  particular,  the  test  (12),  with 


R  _  %  R  =  2._58  (16) 
°n  °  /M 

will  also  give  P(T|T)  =  0.95  and  P(T|0)  =  C.Ol  when  D  =  MDQ,  Since  S/N  <<  1 
when  0  is  near  the  MDQ  level  for  the  N2O  detection  band  used  in  the  field 


tests  (and  tor  all  detection  bands  analyzed  in  Ref.  1),  the  tests  (61  and 
(12)  will  give  the  same  results.* 

The  correlat'cn  coefficient  would  be  a  better  detection  parameter  than 
D'  when  S/N  is  large  and  there  are  spectral  interferences  corresponding  to 
very  large  S/N. 

2 . 2  Calculation  of  P ( T  j  T)  and  P ( T ) 0 ) 

If  D  ^  MDQ.  anJ/or  the  threshold  (R  or  R^)  is  not  set  to  the  value  given 
by  Eq.  (9)  or  (16),  the  detection  and  false  detection  probabilities  can  be 
different  f'om  0.95  and  0.01.  If  detection  is  based  on  the  correlation 
coefficient  p  (and  if  S/N  <<  1),  the  general  equations  are 


P (T |T) 


g(x)  dx 


07) 


oo 

P(T ;  0)  -  2  f  g(x)  dx 

J 

q 


where  g(xl  =  ( 1  /  V2v)  exp(-xV2)  Is  the  unit  normal  density  function  and 


(18) 


[f • ("S/N")  -  R] 

(19) 

R  /M  . 

(20) 

The  ratio  "S/N"  (which  was  denoted  Z  in  Ref.  1)  is  given  by  Eq.  (5).  The 
parameter  f  <_  1  i s  defined  as  the  correlation  coefficient  between  the 
reference  spectrum  and  a  very  high  S,' N  measured  contrast  spectrum.  Generally 
a  value  less  than  unity  indicates  that  the  reference  spectrum  is  not  Perfect 
and/or  that  the  spectrometer  system  is  not  perfectly  calibrated.  Note  that 
P'T|T)  will  be  less  than  0.95  if  f  <  1,  even  if  D  =  MDQ  and  R  is  given  by 
Eq.  ( 1 6)  .  The  integrals  in  (17)  and  (18)  are  tabulated  in  mathematical 
reference  books. 


The  MOQ's  given  in  Ref.  1  are  correct.  Approximate  values  could  have  been 
calculated  from  standard  tables  of  the  normal  distribution,  rather  than  the 
numerical  procedure  defined  in  Ref.  1. 
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If  detection  is 
t i ons  for  P ( T  | '  )  ana 


based  on  the  estimate  D‘  of  D  (llq.  6),  the  above  equa- 
P (T  j  0)  are  valid,  provided  p  and  q  are  redefined: 


P  +  P'  =  (f • |»|  -  Rd]/o0,  , 


(21) 


oQ,  -  /2  NESR/(/H  o  )  , 


(22) 


q  -  q'  =  W  * 


(23) 


2  2  "1 

The  detectable  amounts  D'  In  unit  of  (molecules/cm  ) (W/cm-sr-cm  )  can 
be  converted  to  volume  concentration  C  in  ppmV  from  the  formula 

D  '  10°  Tg*k'  1013  1 

'  SB  2.69  k  I0’9  273 

V 

where  p^  is  the  gas  (ambient)  pressure  and  L  is  the  path  length  through  the 
gas . 
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SECTION  3 

COMPUTER  SIMULATION 


The  equations  for  predicting  system  performance  are  based  on  the  assump¬ 
tion  that  ::he  system  spectral  noise  is  Gaussian  (and  white).  It  w a.,  felt 
that  this  assumption  should  be  tested  before  performir^  the  computer  simula¬ 
tion  and  field  data  analyses. 

USU  >rovided  noise  spectra  obtained  by  differencing  spectral  measure¬ 
ments  of  a  blackbody.  The  FTS  system  observed  the  blackbody  at  very  close 
range  t  minimize  any  effects  of  time-varying  atmospheric  absorption.  ARC 
used  the  standard  chi-square  test  for  goodness-of -f i t  to  test  the  hypothesis 
that  the  spectrum  noise  samples  obey  a  normal  distribution.  The  data  samples 
were  spaced  approximately  0.06  <m  *  apart  and  the  FTS  resolution  (full  width 
at  half  maximum)  was  approximately  0.132  cm  Alternate  samples  were  dis¬ 
carded  in  the  tests,  so  that  the  noise  samples  had  an  effective  spacing  of 
~0.12  cm  '  and  could  be  considered  uncorrelated. 

Initially,  the  entire  noise  spectrum  from  2170  cm  '  io  2260  cm  '  was  used, 
The  first  series  of  tests  involved  21  spectra,  each  consisting  of  7^7  usable 
samples.  The  chi-square  test  with  four  degrees  of  freedom  and  a  significance 
level  of  five  percent  was  applied  to  each  spectrum  individually.  All  21  tests 
foiled,  and  by  a  large  margin. 

Figure  5  shows  why  the  tests  fai:ed:  The  noise  spectrum  obviously  is 
not  white,  and  as  a  result,  the  set  of  noise  samples  is  non-normal.  The 
spectral  noise  is  not  white  because  it  has  been  calibrated;  i . e  ,  the  noise 
voltage  has  been  divided  by  the  spectral  responsivity  of  the  FTS  system.  The 
calibrated  noise  is  large  near  the  edges  of  the  detec'ion  band  because  the 
system  includes  a  cooled  circular  variable  filter  (CVF)  which  has  low  trans¬ 
mission  near  the  band  edges.  The  purpose  of  the  CVF  is,  of  course,  to  reduce 
the  total  photon  flux  on  the  detector;  i.e.,  to  minimize  the  system  NESR. 
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Figure  5-  Example  of  a  (calibrated)  system  noise  spectrum  obtained 
by  differencing  two  measured  spectra  of  a  blackbody. 

The  noise  should  be  white  if  it  is  not  calibrated.  In  fact,  we  realized 
from  these  initial  tests  that  it  would  be  best  to  base  detection  on  the 
analysis  of  uneal ibrated  spectral  data.  The  FTS  spectral  responsivity  func¬ 
tion  should  be  applied  instead  to  the  reference  contrast  spectrum;  i.e.,  the 
reference  spectrum  should  be  multiplied  by  the  spectral  responsivity.  Then 
detection  would  be  based  on  the  correlation  of  two  spectra  with  units  of  volts 
rather  than  spectral  radiance  units.  The  multiplicative  calibration  procedure 
was  used  in  the  analyses  of  field  data  reported  in  Section  4. 

Rather  than  attempt  to  "unca 1 i brate"  the  available  noise  spectra,  we 
narrowed  the  spectral  region  to  2190-2230  cm  ' ,  over  which  the  noise  is 
approximately  white  (see  Fig.  5).  Some  of  the  noise  spectra  had  a  non-zero 
spectrally  varying  average  value,  indicating  that  the  system  spectral  respon¬ 
sivity  had  varied  during  or  between  the  pair  of  blackbody  spectral  measurements. 
Noise  spectra  having  this  behavior  in  the  narrowed  spectral  region  were 
discarded.  Twenty-six  spectra  of  333  samples  each  were  individually  subjected 
to  a  pair  of  chi-square  tests: 

TEST  A:  Two  degrees  of  freedom,  significance  level  five  percent 

TEST  B:  Two  degrees  of  freecom,  significance  level  one  percent 
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It  was  necessary  to  decrease  the  number  or  degrees  of  freedom  from  four  to 
two  because  of  the  reduced  sample  size.  Twenty-three  of  the  2o  spectra 
passed  Test  A.  All  of  fhe  spectra  passed  Test  B.  It  is  concluded  that  i he 
noise  is  normal  (and  white)  over'  the  region  of  uniform  spectral  response,  and 
is  also  normal  (and  white)  over  the  entire  2170-2260  cm  *  detection  hand  if 
the  measurements  are  not  calibrated. 

The  remainder  of  this  section  describes  a  computer  simulation  that  vali 
dated  the  performance  statistics  predicted  by  the  model  described  in  Section  2. 
It  was  recognized  that  while  the  field  tests  might  provide  an  impressive  real- 
world  demonstration  of  the  method,  their  number  would  not  be  sufficient  to 
compile  meaningful  statistics.  it  was  practical  to  simulate  several  thousand 
tests  on  the  computer  whereas  the  planned  number  c  field  tests  was  only  100. 

The  simulation  consisted  of  repeatedly  adding  computer-generated  noise 
to  a  high-S/N  contrast  spectrum  measured  by  'JSU.  This  rerulted  in  an  artifi¬ 
cial  data  base  that  could  be  analyzed  in  the  same  way  as  real  data.  We  first 
analyzed  the  high-S/N  contrast  spectrum  to  estimate  the  detectable  quantity  D 
(from  £q .  2).  Knowing  the  rms  spectral  signal  ~  o^D  to  high  accuracy,  we 
could  adjust  the  rms  value  on  of  the  added  noise  to  obtain  a  particular  value 
of  s  ignal -to-noi  se ,  "S/N"  =  (o^D/cf^)  ^9*  fr°r  var'ous  choices  of  "S/N" 

and  threshold  Fv ,  we  could  forecast  the  detection  and  false  detection  proba¬ 
bilities  (Eqs.  17  through  20).  The  "data"  was  analyzed  (according  to  F.qs.  11 
and  12),  and  the  fraction  of  successful  detections  was  compared  to  the  forecast 
detection  probability.  The  computer-generated  noise  was  analyzed  in  the  same 
way  to  obtain  a  comparison  of  actual  and  forecast  false  detection  probabilities. 

Figure  6  shows  the  high-S/N  measured  contrast  spectrum  on  a  normalized 
(0  to  1.0)  scale.  The  spectrum  is  calibrated,  and  has  been  truncated  to  the 
spectral  region -"2192-2252  cm  ' .  The  non-white  character  of  the  noise  is  not 
important  since  the  S/N  is  reduced  considerably  in  the  simulation.  The  effec¬ 
tive  sample  spacing  is  —0.12  cm  '  and  the  effective  number  of  samples  M  is 
500 .  * 

All  spectrum  plots  in  this  report,  however,  show  computed  or  measured  spectra 
at  —  0.132  cm  '  resolution,  with  points  spaced  at  0.0b  cm  '. 
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Figure  6.  Truncated  USU  high-S/N  contrast  spectrum  used  in 

comouter  simulation  (spectra)  range  ~  2192-2292  cm  ) 

it  was  determined  that  the  high-S/N  measured  contrast  spectrum  was  not 
perfectly  correlated  with  the  reference  spectrum:  the  correlation  coeffi¬ 
cient  was  0.902.  This  value  of  f  was  used  in  Eq.  (19)  to  calculate  the 
predicted  detection  probability. 

The  reference  spectrum  was  computed  using  AFGL's  FASC0D1  computer  code^ 
and  meteorological  data  supplied  by  USU.  The  instrument  spectral  line  shape 
(ILS)  corresponding  to  the  i nterferogram  apodization  used  by  USU  is  called 
the  Happ-Genzel  line  shape. ^  It  was  necessary  to  incorporate  this  particular 
ILS  into  one  of  FASCODl's  subroutines. 

The  results  of  the  simulation  are  shown  in  Table  1.  The  four  cases 
represent  selected  combinations  of  "S/N"  and  threshold  R  which  correspond  to 
a  wide  range  of  predicted  detection  and  false  detection  probabilities.  The 
actual  probabilities  are  based  on  analysis  of  1,000  s igna 1 -pi us-noi se  spectra. 
The  table  gives  two  values  for  each  case  because  the  analysis  was  performed 
twice  for  independent  noise  sets.  The  table  includes  a  comparison  of  pre¬ 
dicted  and  actual  relative  uncertainty  in  the  inferred  detectable  quantity, 

Op  ./D  [The  predicted  value  of  a  /D  is  ( 1  / /m)  ("N/S1')  ,  from  Eqs.  (4)  and  (5)]- 
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As  the  table  shows,  tnere  is  excellent  agreement  between  the  actual  and 
predicted  values,  except  when  the  predicted  P ( T | 0 )  is  much  less  than  one  per¬ 
cent  (Case  2).  There  is  part'cularly  good  agreement  for  the  first  case, 
which  is  closest  to  the  dOQ  performance  level. 

Note  that  we  have  really  tested  the  method  at  a  level  very  close  to  the 
MDQ,  for  if  the  high-S/N  USU  contrast  spectrum  and  reference  spectrum  were 
perfectly  correlated  (f  =  1  rather  than  0,902  in  Eq.  (19)),  the  predicted 
P (T  |  T)  for  the  first  case  could  be  0.9^  rather  than  0,88.  The  actual  P(T|T) 
would  also  be  higher,  and  in  good  agreement  with  the  prediction.  In  the 
actual  practice  of  the  detection  method,  "perfect"  reference  spectra  could  be 
obtained  from  high  S/N  laboratory  measurements  using  the  actual  FTS  detection 
system;  i.e.,  determination  of  the  reference  spectrum  could  be  part  of  the 
system  calibration  procedure. 
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SECTION  4 

ANALYSIS  OF  FIELD  MEASUREMENTS 


The  problem  of  time-varying  spectral  responsivity  referred  to  in  the 
previous  section  was  encountered  again  during  analysis  of  the  low-S/N  field 
measurements.  When  no  target  N^O  is  present  the  contrast  radiance  spectrum 
should  be  just  noise,  but  if  the  spectral  responsivity  has  changed  between 
the  two  measurements,  their  difference  (the  spectral  contrast)  will  contain 
an  additional  component  with  a  slow  spectral  variation.  This  component  will 
result  in  a  false  detection  if  it  is  comparable  to  the  signal  level  and 
if  it  resembles  even  slightly  the  band  contour  of  the  reference  spectrum. 
Similarly,  the  varying  responsivity  can  result  in  missed  detections  and/or 
large  errors  in  the  inferred  detectable  quantity  when  the  target  N^O  is  present 
These  effects  were  observed  in  the  analysis  of  the  first  several  sets  of  field 
data . 


The  remedy  was  simple.  USU  modified  the  FTS  system  to  bypass  the  cooled 
CVF  in  the  optical  train,  and  the  problem  disappeared. 

Estimates  of  the  effective  rms  system  noise  on  =  /l  NESR  were  needed  to 
set  the  detection  threshold,  to  estimate  the  MDQ,  and  to  determine  the  expected 
uncertainty  in  the  inferred  detectable  amounts.  A  single  estimate  was  obtained 
for  each  series  of  runs  by  computing  the  average  of  the  rms  uncalibrated 
spectral  radiance  difference  between  consecutive  pairs  of  background  runs. 

For  each  series,  the  value  of  on  thus  obtained  was  approximately  1.2  mV. 

Figure  7  shows  individual  on  values  determined  from  consecutive  background 
spectrum  pairs  for  five  series  of  runs  made  in  January  and  Februai >  1931. 

The  fact  that  the  variation  in  on  over  a  series  of  runs  slightly  exceeds  the 
95  percent  confidence  limits  on  this  estimate  indicates  that  the  background 
temperature  and/or  FTS  responsivity  may  have  varied  slightly  during  a  series. 
However,  the  variations  are  extremely  small  compared  to  those  observed  when 
the  CVF  was  used  in  the  system. 
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Tp  CONF I DF.NCr  INTERVAL  ON  NOISE  EST1  MATES 


])0l _ I _ L _ I _ 1 _ i _ 1 _ I _ 1 _ J 

1  2  3  4  S  6  7  3  9  10 

RON  NUMBER 

Figure  7.  Noise  estimates  from  differences  of  successive  background  runs 

A  typical  N^O  reference  contrast  spectrum  was  shown  earlier  (Fig.  3). 

For  the  January  and  February  series  of  field  tests  the  standard  deviation 
of  the  computed  reference  spectrum  varied  between  ’.33  x  10  ' ^  and 
l.A  x  10  (cm  /mol ec) (vol ts/W/cm  -sr-cm  ). 

Figure  8  shows  a  typical  FASCODl  computed  transmittance  si  actrum  for  the 
USU  273_meter  path.  The  major  absorptions  are  by  atmospheric  N20  and 
The  transmittance  spectrum  is  included  in  the  "calibrated"  reference  spec¬ 
trum  TvagvRv>  where  is  the  system  spectral  respons i v i ty .  Spectral  elements 
with  x^  <  0.1  were  excluded  from  the  detection  processing.  For  the  January 
and  February  tests  the  effective  number  of  spectral  elements  M  was  between 
709  and  8 1 3 • 

Figure  9  is  a  typical  radiance  spectrum  computed  for  the  USU  273-meter 
path.  The  figure  gives  an  indication  of  how  the  measured  target  and  back¬ 
ground  spectra  would  look  if  they  were  calibrated. 

Tables  A- 1  through  A-5  in  Appendix  A  reproduce  the  computer  printouts  of 
the  data  analysis  results  for  the  January  and  February  field  tests.  These 
results  are  summarized  in  Table  2.  The  second  column  of  the  table  gives  the 
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WAVELENGTH  (MICRON) 
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Figure  8. 


WAVENUMBER  (CH-I) 

Typical  FASCODi  computed  spectral  transmiss 
*or  the  USU  273_meter  path 

WAVELENGTH  (MICRON) 

4.53  4.50 


?2ZCJ 

**  WAVENUMBER  (  CM-  ]  ) 

Figure  9.  Typical  FASCODI  computed  spectral  radiance  for 
the  USU  273-meter  path 

number  of  detections,  and  the  third  the  computed  MDQ.  The  remaining  columns 
give  the  average  of  the  inferred  detectable  quantity  for  the  series  of  12 
runs,  the  standard  deviation  of  these  estimates,  and  the  predicted  standard 
deviation,  all  normalized  to  the  MDQ.  The  results  of  the  19  February  test, 
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in  which  the  flow  rate  was  purposely  varied,  was  already  described  in 
Section  1.  The  second  part  of  the  table  compares  the  average  D'  converted  to 
N.,0  concentration  and  the  measured  concentration  for  two  of  the  field  tests. 

It  is  seen  that  the  agreement  between  inferred  and  measured  concentra¬ 
tions  is  essentially  within  the  calculated  r  one-sigma  uncertainty  in  the 
inferences.  Also,  the  actual  relative  standard  deviations  agree  well  with 
the  estimated  value,  except  for  the  "February"  data  set.  The  small  differences 
(0.3  compared  to  0.2k)  are  probably  due  to  background  varia  ons  or  very  small 
responsivity  changes.  ! 
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SECTION  5 

ESTIMATION  OF  GAS  TEMPERATURE  AND  COLUMN  THICKNESS 


The  temperature  Ts^c  of  the  background,  if  :t  is  assumed  to  radiate  as 
greybody  or  blackbody,  can  be  estimated  by  careful  analysis  of  the  measured 
background  spectrum  and/or  target-plus-background  spectrum.  Depending  on  the 
target  species,  it  may  be  desirable  to  use  a  spectral  region  different  from 
the  detection  band,  i.e.,  an  atmospheric  window  for  which  x  ~  1  in  many  or 
all  spectral  elements.  Once  T^ ^  is  known,  the  target  gas  temperature  T^ 
and  column  thickness  u  can  be  estimated  by  a  more  detailed  analysis  of 
spectra  measured  in  the  detection  band. 


Let  a  denote  AN  +  NOISE  and  b  a  reference  contrast  spectrum.  If  it 

v  v  v  v 

is  desirable  to  subtract  the  average  values  of  these  two  contrast  spectra 
before  correlation  (e.g.,  if  the  background  is  known  to  have  large  spatial  or 
temporal  variations  on  the  scale  of  the  IFOV  separation  and  spectrum  measure¬ 
ment  time),  let  a  and  b  be  defined  as  the  remainders.  Since  T  is  to  be 
v  v  g 

inferred  it  is  appropriate  to  redefine  the  reference  contrast  spectrum  b  to 
include  all  dependences  on  T  : 


i  a„ 

v  gv 


t  a 
v  gv 


(VaVVTsfc> 


The  new  reference  spectrum  b 
a  narticular  guess  (T')  for 
that  minimize  the  rms  error 


v 

T 

g 

i  n 


car  be  computed  (as  before,  via  FASC0D1)  using 
.  The  object  is  to  obtain  estimates  u'  and  [' 
the  equation 


a 

v 


u'bv(T'’Tsfc> 


(25i 


which  is  merely  Eq .  (1)  written  in  different  notation.  The  least  squares 
minimization  procedure  for  determining  u'  and  T'  is  nonlinear,  whereas  the 
one  for  estimating  D  was  linear  (within  the  stated  approximations). 

By  formally  minimizing  the  mean-square  difference  between  the  right- 
an.d  left-hand  sides  of  Eq .  (25)  with  respect  cc  both  u'  and  T',  we  obtain  a 
pair  of  solutions  for  u‘: 
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£avbv(T'-Tsfc>/£bv2(T"’Tsfc.> 


(26a) 


9b 

IX  ^ 


(26b) 


These  two  prescriptions  for  obtaining  u'  can  be  evaluated  for  a  systematically- 
varied  guess  T',  until  they  give  the  same  result.  The  resultant  u'  and  T"  are 
the  least-squares  best  solutions  for  u  and  T  . 

The  solution  method  described  is  quite  straightforward.  However,  anal¬ 
ytical  estimation  of  the  errors  in  the  solutions  due  to  system  no»se  in  a^ 
involves  some  tedious  algebra.  We  performed  a  simplified  analysis  that  assumes 
small  errors.  We  will  present  only  the  results  of  the  analysis. 


The  rms  relative  error  6u/u  in  inferred  column  thickness  u  and  the  rms 
error  fiT^  in  inferred  target  gas  temperature  Tg  are  given  by 


/MZ  ( 


u. 


a  1  nAB 

3T9 


where 


9  /MZ  (2k) 


uob/(/2  NESR)  , 

I  in,  1  - 

AT->0  /*t\2 


The  derivative  in  Eq .  (27)  is  9 1 nAB ^ ( T , T ^ ) /3T  evaluated  at  T  *  T  .  The  rms 
value  of  b^  (new  definition)  is  denoted  o^;  hence  Z  is  essentially  the  same 
as  the  "S/N",  given  by  Eq.  (5).  The  parameter  k  may  be  defined  as  the 
coefficient  of  the  quadratic  term  in  the  Taylor  expansion  of  1  —  p ,  where  p  is 
the  correlation  coefficient  of  the  reference  spectra  b^(T^)  and  b^(Tq+AT). 
Essentially,  k  is  a  measure  of  the  sensitivity  of  the  spectral  contour  (rota¬ 
tional  line  intensity  distribution)  of  the  target  gas  absorption  band  to  a 
change  in  the  gas  temperature.  Note  that  6T^  is  small  for  large  k. 
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It  should  be  mentioned  that  the  random  variates  u'  and  T',  whose 
standard  deviations  are  given  by  Eqs.  (27)  and  (28),  are  partially  correlated. 
That  is,  the  uncertainty  or  "noise"  in  u'  is  the  sum  of  two  partially  corre¬ 
lated  normal  variates,  one  of  which  is  proportional  to  the  "noise"  in  T'. 

Equation  (27)  is  always  an  overestimate  of  the  rms  error  6u/u  since  it  assumes 
perfect  positive  correlation  between  the  components  of  u'. 

Equations  (27)  and  (28)  were  evaluated  for  conditions  typical  of  the  US IJ 
field  tests.  The  values  T^  =  295K  and  T^^c  =  309 K  were  used,  and  it  was 

assumed  that  K'„0  is  present  at  the  MDQ  level.  It  was  also  assumed  that  T  is 

L  ft;V  9 

the  same  as  the  atmospheric  temperature.  The  calculations  gave 


Su/u  ~  4.0 

6T  cs  6>3K 
9 


(3D 


which  are  absurdly  large  errors.  Also,  they  are  probably  not  accurate  error 

estimates  since  Eqs.  (27)  and  (28)  were  derived  on  the  assumption  of  small 

errors.  They  show,  however,  that  the  errors  would  be  of  the  order  of  <Su/u  2*  0 

and  ST  ~  6.3K  if  N~Q  were  present  at  ten  times  the  MDQ  level, 
y  ^ 

These  results  confirmed  our  expectation  that  u  and  Tg  cannot  be  deter¬ 
mined  accurately  from  the  target  gas  contrast  spectrum  alone  when  the  gas  is 
present  at  the  MDQ  level.  If  some  other  "tracer"  gas  were  present  in  such 
large  quantity  that  T^  could  be  determined  almost  exactly  from  its  spectrum, 
then  the  error  6u/u  would  be 


5u  _ ^  1 

U  /MZ 

which  has  the  value  0.24  for  N^Q  at  the  MDQ  level. 


(32) 


Knowledge  of  the  covariance  matrix  would  be  useful  in  defining  a  two-dimers io 
decision  region  that  would  maximize  detection  probability  and  minimize  false 
detection  probability.  It  would  also  allow  accurate  quantification  of  the 
errors  6u  and  5T  We  did  not  evaluate  the  correlation  matrix  because  of  the 
great  amount  of  Tabor  and  computing  time  it  would  have  required. 

And  that  it  undergoes  the  same  variation  as  the  atmospheric  temperature. 

This  assumption  effects  the  value  of  k  calculated  from  Eq.  (30). 

The  error  in  Sv  is  also  strongly  dependent  on  the  temperature  difference 
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